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Abstract 

The status of the evaluation of the MSSM Higgs sector is reviewed. The phe- 
nomenological impact of recently obtained corrections is discussed. In particular it 
is shown that the upper bound on nih within the MSSM is shifted upwards. Conse- 
quently, lower limits on tan /5 obtained by confronting the upper bound as function 
of tan P with the lower bound on nih from Higgs searches are significantly weakened. 
Furthermore, the region in the M^-tan /3-plane where the coupling of the lightest Higgs 
boson to down-type fermions is suppressed is modified. The presently not calculated 
higher-order corrections to the Higgs-boson mass matrix are estimated to shift the 
mass of the lightest Higgs boson by up to 3 GeV. 
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1 Introduction 



A crucial prediction of the Minimal Supersymmetric Standard Model (MSSM) [Ij is the 
existence of at least one light Higgs boson. The search for this particle is one of the main 
goals at the present and the next generation of colliders. Direct searches at LEP have already 
ruled out a considerable fraction of the MSSM parameter space, and the forthcoming high- 
energy experiments at the Tevatron, the LHC and a future Linear Collider (LC) will either 
discover a light Higgs boson or rule out the MSSM as a viable theory for physics at the weak 
scale. Furthermore, if one or more Higgs bosons are discovered, their masses and couplings 
will be determined with high accuracy at a future LC. Thus, a precise knowledge of the 
dependence of the masses and mixing angles of the MSSM Higgs sector on the relevant 
supersymmetric parameters is of utmost importance to reliably compare the predictions of 
the MSSM with the (present and future) experimental results. 

We recall that the Higgs sector of the MSSM |2| consists of two neutral CP-even Higgs 
bosons, h and H {rrih < mfj), the CP-odd A boson, and two charged Higgs bosons, H"^. At 
the tree-level, mh,tree and mH,tTee can be calculated in terms of the Standard Model (SM) 
gauge couplings and two additional MSSM parameters, conventionally chosen as Ma and 
tan/?, the ratio of the two vacuum expectation values (tan/5 = V2/V1). The two masses are 
obtained by rotating the neutral CP-even Higgs boson mass matrix with an angle a, 

2,tree ^ ( siu^ /? + M| COS^ /3 - (M^ + M|) siu /? COS /? \ 

^^ss^ 1^ -(Ml + M|) sin ^ COS /? Ml cos^ /3 + M| sin^ /5 )' ^ ' 



with a satisfying 



tan2a = tan2/? ^^_^f , -- < a < 0. (2) 



In the Feynman diagrammatic (FD) approach the higher-order corrected Higgs boson 
masses are derived by finding the poles of the h, if-propagator matrix whose inverse is given 

by 

-1 _ J - ^H, tree + ^Hh{P^) ^hH{p^^ ^ 



^ ^""^ I tHH{p') p'-<,.ee + ^HH{p'))' ^ ^ 

where the S(p^) denote the renormalized Higgs-boson self-energies, p being the momen- 
tum flowing on the external legs. Determining the poles of the matrix Aniggs in eq. (jH)) is 
equivalent to solving the equation 

P^ - ^h,tree + Sm(/)] " m'^H,tree + ^Hnip'^)] " [Sft/f(p^)]^ = . (4) 

The status of the available results for the self-energy contributions to eq. © can be 
summarized as follows. For the one-loop part, the complete result within the MSSM is 
known The by far dominant one-loop contribution is the 0{at) term due to top and 

stop loops (at = hf/{4:TT), ht being the superpotential top coupling). Concerning the two-loop 
effects, their computation is quite advanced and it has now reached a stage such that all the 
presumably dominant contributions are known. They include the strong corrections, usually 
indicated as 0{atas)i and Yukawa corrections, 0{a1)i to the dominant one-loop 0{at) term, 
as well as the strong corrections to the bottom/sbottom one- loop 0{ab) term (a^ = /i^/(47r)), 
i.e. the 0{ahas) contribution. The latter can be relevant for large values of tan /3 . Presently, 
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the 0{atas) [313131101111], 0{af) and the 0{abas) [H contributions to the self- 

energies are known for vanishing external momenta. In the (s) bottom corrections the all- 
order resummation of the tan (3 -enhanced terms, 0{ab{as tan is also performed [T^IITB] . 

In this paper we are going to present, in view of the recent achievements obtained in 
the knowledge of the two-loop corrections, updated results for various quantities of physical 
interest. In our analysis we employ the latest version of the Fortran code FeynHiggs fr7\\lS\ . 
namely FeynHiggsl.3, that evaluates the MSSM neutral CP-even Higgs sector masses and 
the mixing angle. 

The paper is organized as follows: in Section |21 we give a brief summary of recent theoret- 
ical improvements in the MSSM Higgs sector and describe the corresponding modifications 
in the code FeynHiggs in order to comprise all existing higher-order results. We compare 
the results for rrih obtained employing different approximations in the treatment of the two- 
loop corrections. In Section |21 we present the results for the phenomenological implications 
of our improved knowledge of the two-loop contributions to the Higgs boson self-energies. 
Section m contains a discussion of the presently unknown contributions in the prediction for 
the CP-even Higgs-boson masses of the MSSM, and an estimate of their possible numerical 
importance is given. Finally, in Section we draw our conclusions. 

2 Recent improvements in the MSSM Higgs sector and 
implementation in FeynHiggs 

In order to discuss the impact of recent improvements in the MSSM Higgs sector we will 
make use of the program FeynHiggs |17| ll8j. which is a Fortran code for the evaluation 
of the neutral CP-even Higgs sector of the MSSM including higher-order corrections to the 
renormalized Higgs boson self-energies. The original release included the well known full one- 
loop corrections jnUHElEI, the two-loop leading, momentum-independent, 0{atas) correction 
in the top/stop sector [Hlinilll] as well as the two- loop leading logarithmic corrections at 
Ola'f) [T^I2Uj . By two-loop momentum-independent corrections here and hereafter we mean 
the two-loop contributions to Higgs boson self-energies evaluated at zero external momenta. 
At the one-loop level, the momentum-independent contributions are the dominant part of 
the self-energy corrections, that, in principle, should be evaluated at external momenta equal 
to the poles of the h, if-propagator matrix, eq. 

The new version of FeynHiggs contains a modification in the one-loop part due to a 
different renormalization prescription employed and includes also several new corrections to 
the Higgs boson self-energies that have recently been calculated. These changes are described 
in the following subsections. 

With the implementation of the latest results obtained in the MSSM Higgs sector, Feyn- 
Higgs comprises all available higher-order corrections and thus allows the presently most pre- 
cise prediction of the masses of the CP-even Higgs bosons and the corresponding mixing an- 
gle. The latest version of FeynHiggs, FeynHiggsl.3, can be obtained from www . f eynhiggs . de. 

2.1 Hybrid renormalization scheme at the one-loop level 

FeynHiggs is based on the FD approach with on-shell renormalization conditions ^9j. This 
means in particular that all the masses in the FD result are the physical ones, i.e. they 
correspond to physical observables. Since eq. (jH) is solved iteratively, the result for nih 
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and rriH contains a dependence on the field- renormalization constants of h and H, wliicli 
is formally of higher order. Accordingly, there is some freedom in choosing appropriate 
renormalization conditions for fixing the field-renormalization constants (this can also be 
interpreted as affecting the renormalization of tan/5). Different renormalization conditions 
have been considered in the literature, e.g. (S' denotes the derivative with respect to the 
external momentum squared): 

1. on-shell renormalization for S^, S^, S^, T^az, and 6vi/vi = 6V2/V2 [HI 

2. on-shell renormalization for S^, ^a, ^az, and 6vi = 5t>j,div, i = 1,2 0121] 

3. on-shell renormalization for S^, , DR renormalization (employing dimensional re- 
duction |22]) for SZhySZn and tan/? jlHj. 

Versions of FeynHiggs previous than the 1.2 release were based on type-1 renormalization 
conditions, thus requiring the derivative of the A boson self-energy. The latest versions 
employ instead the hybrid DR/on-shell type 3 conditions ^Hj- The choice of a DR defini- 
tion^ for SZh^SZn and tan/? requires to specify a renormalization scale at which these 
parameters are defined. This scale is set to in FeynHiggs, but can be changed by the 
user. These new renormalization conditions lead to a more stable behavior around thresh- 
olds, e.g. Ma = ^rrit, and avoid unphysically large contributions in certain regions of the 
MSSM parameter space (a detailed discussion can be found in Ref. p3j; see also Ref. [21! )• 

2.2 Two-loop 0{a^) corrections 

In the previous version of FeynHiggs the two-loop 0{af) corrections were implemented in 
a simplified form, taken over from the result of P^ I2Uj. obtained with the renormalization 
group (RG) method. The latter result included the two-loop leading logarithmic corrections, 
proportional to the square of t = log {M^^o^y/mf) (where Msusy is a common scale for the 
soft SUSY-breaking parameters), while the next-to-leading logarithmic terms, linear in t, 
were not totally accounted for. 

Recently, the two-loop 0{a^) corrections in the limit of zero external momentum became 
available, first only for the lightest eigenvalue, rrih, and in the limit Ma 3> Mz ^2]; then for 
all the entries of the Higgs propagator matrix for arbitrary values of Ma They were 
obtained in the effective-potential approach, that allows to construct the Higgs boson self- 
energies, at zero external momenta, by taking the relevant derivatives of the field-dependent 
potential. In this procedure it is important, in order to make contact with the physical 
Ma, to compute the effective potential as a function of both CP-even and CP-odd fields, 
as emphasized in Ref. In the evaluation of the 0{af) corrections, the specification 

of a renormalization prescription for the Higgs mixing parameter /i is also required and 
it has been chosen as DR. In FeynHiggsl.3, which includes the complete two-loop 0{af) 
corrections, the corresponding renormalization scale is fixed to be the same as for 6Zh, SZh 
and tan /5. 

The availability of the complete result for the momentum-independent part of the 0{af) 
corrections allows us to judge the quality of results that incorporate only the logarithmic 
contributions. In Fig. Q we plot the two- loop corrected rrih as a function of the stop mixing 

"'^A minimal definition for tan/3 and the field renormalization constants was already employed in Ref. 0], 
in the analysis of the one-loop (s)top-(s)bottom contribution to the Higgs boson self-energies. 
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Figure 1: Two- loop corrected mh as a function of Xf in various steps of approximation. 
The relevant MSSM parameters are chosen as tan/? = 3 , m^^'' = 174.3 GeV, M^^ = M-^^ = 
Ma = /i = 1 TeV and rrig = 800 GeV. The meaning of the different curves is explained in 
the text. 

parameter, Xt = At—fi/ tan/3, where At denotes the trilinear soft SUSY breaking Higgs-stop 
coupling. We use the following convention for the stop mass matrix, 

A^2 ^ / + + COS 2f3 (i - |s2^)M| mtXt 

* mtXt Mf^+mf + I cos2(3 slrM^ 

For simplicity, the soft SUSY breaking parameters in the diagonal entries of the stop mass 
matrix, Mj^, M^^, are chosen to be equal, M^^ = M^^ = Msusy- For the numerical analysis 
MsusY as well as Ma and fi, are chosen to be all equal to 1 TeV, while the gluino mass 
is rrig = 800 GeV, and tan/? = 3. If not otherwise stated, in all plots below we choose 
the trilinear couplings in the stop and sbottom sectors equal to each other, Ab = At, and 
set M2 = Msusy, where M2 is the SU(2) gaugino mass parameter (the U(l) gaugino mass 
parameter is obtained via the GUT relation. Mi = 5/3 Sy^/c^^r M2). 

The solid and dotted lines in Fig. Q are computed with and without the inclusion of 
the full 0{at) corrections, while the dashed and dot-dashed ones are obtained including 
only the logarithmic contributions. In particular, the dashed curve corresponds to the result 
of [ini, obtained with the one-loop renormalization group method. The dot-dashed curve. 
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instead, corresponds to the leading and next-to- leading logarithmic terms that, within this 
simplifying choice of the MSSM parameters, can be easily singled out from the complete 
0{af) result. As mentioned above, the two approximate results agree with each other in 
the terms proportional to t^, whereas they disagree in the terms linear in t, giving rise to 
the difference between the dashed and dot-dashed curves. In fact, as long as one is only 
concerned about the leading logarithmic terms, it is correct to use in the renormalization 
group equations (RGE) for the Higgs quartic couplings the one-loop /3-function as is done 
in Ref. pl9|. Moreover, the exact definition of the various parameters in the one- loop part 
is not important. Instead, when the next-to- leading logarithmic terms are examined the 
two-loop /3-function has to be employed, giving rise to additional single logarithmic terms 
at the two-loop order. Furthermore, the results from |Tni are valid under the assumption 
that the one-loop part of the corrections is written in terms of running SM parameters in 
the MS renormalization scheme, whereas the one-loop computation in FeynHiggs employs 
on-shell parameters. As discussed in |20[TU[|23I2S1 , this amounts to a shift in the two- loop 
corrections that also affects the next-to-leading logarithmic part of the 0{af) results. 

From Fig. ^ it can also be seen that, for small Xt, the full 0{af) result is very well 
reproduced by the logarithmic approximation, once the next-to-leading terms are correctly 
taken into account. On the other hand, when Xt is large there are significant differences, 
amounting to several GeV, between the logarithmic approximation and the full result. Such 
differences are due to non-logarithmic terms that scale like powers of Xt/MsvsY- It should 
be noted that for more general choices of the MSSM parameters the renormalization group 
method becomes rather involved (see e.g. ^7\, where the case of a large splitting between 
Mf^ and M^^ is discussed), and a suitable next-to- leading logarithmic approximation to the 
full result is much more difficult to devise. 

2.3 Two-loop sbottom corrections 

Due to the smallness of the bottom mass, the 0{ab) one-loop corrections to the Higgs boson 
self-energies can be numerically non-negligible only for tan /5 3> 1 and sizable values of the 
/i parameter. In fact, at the classical level hb/ht = {nib/mt) tan/3, thus tan/5 3> 1 is needed 
in order to have ab ~ a* in spite of rrib nit. In contrast to the 0{at) corrections where 
both top and stop loops give sizable contributions, in the case of the 0{ab) corrections the 
numerically dominant contributions come from sbottom loops: those coming from bottom 
loops are always suppressed by the small value of the bottom mass. A sizable value of /i is 
then required to have sizable sbottom-Higgs scalar interactions in the large-tan/? limit. 

The relation between the bottom-quark mass and the Yukawa coupling hb, which controls 
also the interaction between the Higgs fields and the sbottom squarks, reads at lowest order 
nib = hbVi/\/2. This relation is affected at one-loop order by large radiative corrections 
jl5j . proportional to hbV2, giving rise in general to tan /^-enhanced contributions. These 
terms proportional to f2, often indicated as threshold corrections to the bottom mass, are 
generated either by gluino-sbottom one-loop diagrams, resulting in O^abOis) corrections to 
the Higgs masses, or by chargino-stop ones, giving 0{abat) corrections. Because the tan/?- 
enhanced contributions can be numerically relevant, an accurate determination of hb from 
the experimental value of the bottom mass requires a resummation of such effects to all 
orders in the perturbative expansion, as described in Ref. 

Concerning the sbottom corrections to the renormalized Higgs boson self-energies, the 
version 1.2.2 of FeynHiggs 'm.chxded the full one-loop contribution, improved by the resumma- 
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Figure 2: The result for the hghtest CP-even Higgs-boson mass in the MSSM, ruh, as 
obtained with the program FeynHiggs is shown as a function of tan (3 for = 120 GeV, 
/i = -1 TeV, = = = m-g = I TeV, At = A^, = 2 TeV. The meaning of the 
different curves is explained in the text. 

tion of the tan /3-enhanced effects in the relation between h}, and according to the effective 
Lagrangian formahsm developed in Ref. [16j. This takes into account the 0{ai,{as taxi [3Y) 
corrections to all orders in n. Numerically this is by far the dominant part of the contribu- 
tions from the sbottom sector. 

Very recently, the complete two-loop, momentum- independent, 0{ahCis) corrections (which 
are not included in the tan/5)") resummation) have been computed [Tl|. The result 

obtained makes use of an appropriate choice of renormalization conditions on the relevant 
parameters that allows to disentangle the genuine two-loop effects from the large thresh- 
old corrections to the bottom mass, and also ensures the decoupling of heavy gluinos. The 
complete momentum-independent 0{ahOis) corrections have now also been implemented in 
FeynHiggsl.3. 

To appreciate the importance of the various sbottom contributions, we plot in Fig. |21the 
light Higgs mass rrih as a function of tan f3. The SM running bottom mass computed at the 
top mass scale, mh{mt) = 2.74 GeV, is used in order to account for the universal large QCD 
corrections. The relevant MSSM parameters are chosen as Ma = 120 GeV, fi = —1 TeV, 
Mf^ ^ Mi^ = Ml =m-g = l TeV, At = Ab = 2 TeV (where is a soft SUSY breaking 



6 



parameter in the sbottom sector). The dot-dashed curve in Fig. |21 includes the full one- loop 
contribution as well as the two-loop 0{atas + af) corrections (the latter being approximately 
tan /^-independent when tan/3 is large). The dashed curve includes also the resummation of 
the tan /^-enhanced threshold effects in the relation between hb and mb- Finally, the solid 
curve includes in addition the complete 0{abas) two-loop corrections of fHl- In the last 
two curves, the steep dependence of rrih on tan/5 when the latter is large is driven by the 
sbottom contributions. We see that, although the tan /^-enhanced threshold effects account 
for the bulk of the sbottom contributions beyond one-loop, the genuine 0{ab(ys) two-loop 
corrections can still shift rrih by several GeV for large values of tan/3 and /z.^ 

3 Phenomenological implications 

The improved knowledge of the two-loop contributions to the Higgs-boson self-energies re- 
sults in a very precise prediction for the Higgs-boson masses and mixing angle with interesting 
implications for MSSM parameter space analyses. In this section we are going to discuss 
possible implications on the upper limit on within the MSSM and on the correspond- 
ing limit on tan/3 arising from confronting the upper bound on with the lower limit 
from Higgs searches. We furthermore investigate the modifications that are induced in the 
couplings of the lightest Higgs boson to the down-type fermions. 

3.1 Limits on ruh and tan/3 

The theoretical upper bound on the lightest Higgs-boson mass as a function of tan jS can be 
combined with the results from direct searches at LEP to constrain tan /?. The diagonaliza- 
tion of the tree- level mass matrix, eq. (Q), yields a value for m^^tree that is maximal when 
M4 3> Mz, in which case ^«^tree — cos^ 2/? , which vanishes for tan/3 = 1. Radiative 
corrections significantly increase the light Higgs-boson mass compared to its tree-level value, 
but still rrifi is minimized for values of tan /? around one. Thus, in principle, the region of low 
tan/? can be probed experimentally via the search for the lightest MSSM Higgs boson j2H]- 
If the remaining MSSM parameters are tuned in such a way to obtain the maximal value of 
ruh as a function of tan j3 (for reasonable values of Msusy and taking into account the ex- 
perimental uncertainties of and the other SM input parameters as well as the theoretical 
uncertainties from unknown higher-order corrections), the experimental lower bound on ruh 
can be used to obtain exclusion limits for tan /?. While in general a detailed investigation of 
a variety of different possible production and decay modes is necessary in order to determine 
whether a particular point of the MSSM parameter space can be excluded via the Higgs 
searches or not, the situation simplifies considerably in the region of small tan/? values. In 
this parameter region the lightest CP-even Higgs boson of the MSSM couples to the Z bo- 
son with SM-like strength, and its decay into a bb pair is not significantly suppressed. Thus, 
within good approximation, constraints on tan j3 can be obtained in this parameter region by 
confronting the exclusion bound on the SM Higgs boson with the upper limit on rrih within 
the MSSM. We use this approach below in order to discuss the implications of the new rrih 
evaluation on tan /3 exclusion bounds. 

^In the region of large tan p the bottom Yukawa coupling receives large corrections, see eq. ((HJl below. 
For the parameter range shown in Fig. |21the quantity A;,, defined in eq. is negative but does not get 
close to —1, |Af,| ^ 0.6. Thus, the perturbative treatment in Fig.|5]seems to be justified. 
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Concerning the upper bound on rrih within the MSSM, the one-loop corrections contribute 
positively to m\. The two-loop effects of 0{atas) and C(a^), on the other hand, enter 
with competing signs, the former reducing m\ while the latter giving a (smaller) positive 
contribution. The actual bound that can be derived depends sensitively on the precise value 
of the top-quark mass, because the dominant one- loop contribution to m\, as well as the 
two-loop 0{atas) term, scale as mf. Furthermore, a large top mass amplifies the relative 
importance of the two-loop 0{af) correction, because of the additional factor. 

In order to discuss restrictions on the MSSM parameter space it has become customary 
in the recent years to refer to so-called benchmark scenarios of MSSM parameters ^29^ .30]. 
The m™*^^ benchmark scenario [211 been designed such that for fixed values of rrit and 
-^susY the predicted value of the lightest CP-even Higgs-boson mass is maximized for each 
value of Ma and tan f3. The value of the top-quark mass is fixed to its experimental central 
value, rrit = 174.3 GeV, while the SUSY parameters are taken as (referring to their on-shell 
values according to the FD result as implemented in FeynHiggs): 

MsusY = 1 TeV, = -200 GeV, Ms = 200 GeV, 

X™ = 2 MsusY, Ah = A, nig = 0.8 MgusY , (6) 

and we will investigate below the case M4 = 1 TeV. 

In Fig. El we plot nih as a function of tan/5 in the m™^^ scenario. The dashed and solid 
curves correspond to the result obtained with the previous (used for the LEP evaluations 
so far |2H1) and the latest (which will be used for the final LEP evaluations ^Tl) versions 
of FeynHiggs, respectively. The two versions, FeynHiggsl.O and FeynHiggsl.3, differ by the 
recent improvements obtained in the MSSM Higgs sector which are described in Sect. |21 
For comparison, also the result obtained with a renormalizat ion-group improved effective 
potential method is indicated. The dotted curve in Fig. |21 corresponds to the code subh- 
poledm 19j,2Q»ilM the m™^^ scenario, for a t mixing parameter ^ = V6Msusy [19,2^. 
It deviates from the result of FeynHiggsl.O by typically not more than 1 GeV for tan/? > 1. 
The LEP exclusion bound for the mass of a SM-like Higgs 123, Mf^ > 114.4 GeV, is shown 
in the figure as a vertical long-dashed line. As can be seen from the figure, the improvements 
on the theoretical prediction described in Sect. O in particular the inclusion of the complete 
momentum-independent 0(a( ) corrections into FeynHiggs, gives rise to a significant increase 
in the upper bound on nih as a function of tan (3. Comparison of this prediction with the 
exclusion bound on a SM-like Higgs shows that the lower limit on tan/? is considerably 
weakened. 

Concerning the interpretation of the results shown in Fig. ^ it should be kept in mind 
that in the m™^^ benchmark scenario rrit and MgusY are kept fixed, and no theoretical 
uncertainties from unknown higher-order corrections are taken into account. In order to 
arrive at a more general exclusion bound on tan/5 that is not restricted to a particular 
benchmark scenario, the impact of the parametric and higher-order uncertainties in the 
prediction for has to be considered j33]. In order to demonstrate in particular the 
dependence of the tan /? exclusion bound on the chosen value of the top pole mass, the dot- 
dashed curve in Fig. |21 shows the result obtained with FeynHiggsl.3 where the top-quark 
mass has been increased by one standard deviation, =5.1 GeV, to rrit = 179.4 GeV, 
and MsusY has been changed from 1 TeV to 2 TeV. It can be seen that in this more general 
scenario no lower limit on tan j3 from the LEP Higgs searches can be obtained. 

Constraints from the Higgs searches at LEP do of course play an important role in regions 
of the MSSM parameter space where the parameters are such that rrih does not reach its 
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Figure 3: The result for the hghtest CP-even Higgs-boson mass, m/i, as a function of tan/3 
in the m™'^'^ scenario. The dotted curve has been obtained with a renormahzation-group 
improved effective potential method, using the code subhpoledm. The dashed curve corre- 
sponds to the result obtained with FeynHiggsl.O, while the full curve shows the result of 
FeynHiggsl.3, where the improvements described in Sect. El are included. The dot-dashed 
curve, also obtained with FeynHiggsl.3, employs rrit = 179.4 GeV and Msusy = 2 TeV. The 
vertical long-dashed line corresponds to the LEP exclusion bound for the SM Higgs boson 
of 114.4 GeV. 

maximum value. Also in this case, however, the remaining theoretical uncertainties from 
unknown higher-order corrections (see Sect. |3] below) have to be taken into account in order 
to obtain conservative exclusion limits. 



3.2 Higgs couplings to fermions 

The tree-level couplings of the lightest CP-even Higgs boson to the up-type and down-type 
SM fermions read, respectively: 

ierriu -cosa _ ierud sina 

~ 2 swMw sin (3 ' ~ 2 swMw cos (3' ^ ' 

where the factors involving a and (3 reflect the changes in the MSSM compared to the SM 
couplings. In the limit of Ma ^ Mz, a ~ (3 — tt/2 so that the SM limit is recovered. 
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For the lightest CP-even Higgs boson the decay to two bottom quarks is usually the main 
decay mode, while the decay to r leptons has usually the second largest branching ratio. 
However, these channels are also most significantly affected by loop corrections, which can 
change the situation described above quantitatively and even qualitatively. The two main 
sources (besides the SM QCD corrections) are the Higgs-boson propagator corrections and 
the corrections modifying the relation between the bottom quark or r lepton mass and the 
corresponding Yukawa couplings. 

Concerning the former, it has been shown analytically [^|33] that the momentum- 
independent contributions coming from the Higgs-boson propagator corrections can be in- 
corporated by replacing the tree- level angle a in eq. ((Tj) with an effective angle ctefr, which 
diagonalizes the Higgs-boson mass matrix including the self-energy corrections evaluated 
at zero external momentum. Due to the effect of the higher-order corrections, Oefr ~ is 
possible, i.e. the coupling of the lightest Higgs boson to the down-type SM fermions can 
vanish. 

The other potentially large corrections to the couplings come from the tan /3-enhanced 
threshold effects in the relation between the down-type fermion mass and the corresponding 
Yukawa coupling [1^], already mentioned in Sect. 12.31 A simple way to take into account 
these effects is to employ the effective lagrangian formalism of Refs. fOlEni! where the cou- 
pling of the lightest Higgs boson to down-type fermions (expressed through the fermion 
mass) is modified according to 



I e rrid sm acs 



eff 2swMw cos/3 



1 



A. 



1 + 



(1 + cot cot P) 



where A^^ contains the tan /3-enhanced terms, and other subleading (i.e. non tan /3-enhanced) 
corrections have been omitted. In the case of the coupling to the bottom quarks, the leading 
contributions to A;, are of 0{as) and 0{at), coming from diagrams with sbottom-gluino and 
stop-chargino loops, respectively. In the case of the r leptons, the leading contributions are 
of 0{ar), coming from sneutrino-chargino loops. The terms containing A^ in eq. (jHl) may 
be relevant for large tan/? and moderate values of Ma- When Ma is much bigger than Mz, 
the product cotaeg cot/5 tends to —1, and the SM limit is again recovered. 

The effects of the higher-order corrections to the couplings of the lightest Higgs boson 
to the down- type fermions appear very pronounced in the "small agff" scenario jSH], corre- 
sponding to the following choice of MSSM parameters: 



mt = 174.3 GeV, Mj 
/i = 2.5MsusY, M2 



SUSY 



800 GeV, 



FD 



500 GeV, Ah 



- K 



m, 



-1100 GeV 

, = 500 GeV . (9) 



In the upper row of Fig. HI we show the ratio sin^aes/ cos^/3 in the M^-tan/? plane, 
evaluated in the "small Oefr" scenario. As explained above, replacing a by aes in the tree-level 
couplings of eq. ((Tj) gives rise to regions in the M^-tan P plane where the effective coupling 
of the lightest Higgs boson to the down-type fermions is significantly suppressed with respect 
to the Standard Model. The region of significant suppression of sin^ aes / cos^ P as evaluated 
with FeynHiggsl.O, i.e. without the inclusion of the sbottom corrections beyond one-loop 
order, is shown in the upper left plot. The upper right plot shows the corresponding result 
as evaluated with FeynHiggsl.3, where the 0{abas) and O {ab{as tan P)^) corrections as well 
as the new 0{a^) ones are included. The new corrections are seen to have a drastic impact on 
the region where sin^ «eff / cos^ P is small. While without the new corrections a suppression 
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Figure 4: Regions of significant suppression of the coupling of h to down-type fermions 
in the M^-tan /5-plane within the "small ctefr" benchmark scenario. The upper left plot 
shows the ratio sin^aefr/ cos^/5 as evaluated with FeynHiggsl.O (i.e. without O^abCis) and 
0(ab(Q;s tan/?)") corrections), while the upper right plot shows the same quantity as evalu- 
ated with FeynHiggsl.3 (i.e. including these corrections). The plot in the lower row shows 
the ratio T{h fc&)MSSM / r(^ bb)sM, i-e. the partial width for h ^ bb normalized to its 
SM value, where the other term in eq. (jH)) and further genuine loop corrections are taken 
into account (see text). 



of 70% or more occurs only in a small area of the M^-tan /3-plane for 20 ^ tan /? < 40 and 
100 GeV < Ma ^ 200 GeV, the region where sin^agfj/ cos^/? is very small becomes much 
larger once these corrections are included. It now reaches from tan /5 ^ 15 to tan /3 > 50, and 
from Ma ^ 100 GeV to Ma ^ 350 GeV. The main reason for the change is that the one- 
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loop 0{ab) corrections to the Higgs-boson mass matrix, which for large tan/3 would prevent 
OefT from going to zero, are heavily suppressed by the resummation of the 0(Q;f,(as tan/5)") 
corrections in the bottom Yukawa coupling. This kind of suppression depends strongly on 
the chosen MSSM parameters, and especially on the sign of /x. 

In order to interpret the physical impact of the effective coupling shown in the upper 
row of Fig. 01 the terms in eq. (jHl) as well as further genuine loop corrections occuring in 
the h —>■ bb process have to be taken into account. The effect of these contributions can be 
seen from the plot in the lower row of Fig. |1[ where the ratio T{h —>■ fefe)MSSM / bb)sM 
is shown, which has been evaluated by including all terms of eq. (jHl) as well as all further 
corrections described in Ref. [331. The region where the partial width r{h bb) within the 
MSSM is suppressed compared to its SM value is seen to be somewhat reduced and shifted 
towards smaller values of Ma as compared to the region where sin^ aes / cos^ /3 is small. 

4 Estimate of the theoretical accuracy of the Higgs- 
boson mass determination 

The prediction for ruh in the MSSM is affected by two kinds of uncertainties, parametric 
uncertainties from the experimental errors of the input parameters and uncertainties from 
unknown higher-order corrections. While currently the parametric uncertainties dominate 
over those from unknown higher-order corrections, as the present experimental error of the 
top-quark mass of about ±5 GeV induces an uncertainty of Am/j ^ ±5 GeV [33], at the next 
generation of colliders rrit will be measured with a much higher precision, reaching the level 
of about 0.1 GeV at an e^e~ LC |2nj- Thus, a precise measurement of rrih will provide a high 
sensitivity to SUSY loop effects and will in this way allow a stringent test of the MSSM, 
provided that the uncertainties from unknown higher-order corrections are sufficiently well 
under control. 

Given our present knowledge of the two-loop contributions to the Higgs-boson self- 
energies, the theoretical accuracy reached in the prediction for the CP-even Higgs-boson 
masses is quite advanced. However, obtaining a complete two-loop result for the Higgs- 
boson masses and mixing angle requires additional contributions that are not yet available. 
In this section we discuss the possible effect of the missing two-loop corrections, and we 
estimate the size of the higher-order (i.e. three-loop) contributions. 

4.1 Missing two- loop corrections 

It is customary to separate the corrections to the Higgs boson self-energies into two parts: 

i) the momentum- independent part, namely the contributions to the self-energies evaluated 
at zero external momenta, which can also be computed in the effective potential approach; 

ii) the momentum-dependent corrections, i.e. the effects induced by the dependence on the 
external momenta of the self-energies, that are required to determine the poles of the h, H- 
propagator matrix. 

All the presently known two-loop contributions are computed at zero external momentum, 
and moreover they are obtained in the so-called gaugeless limit, namely by switching off the 
electroweak gauge interactions. The two approximations are in fact related, since the leading 
Yukawa corrections are obtained by neglecting both the momentum dependence and the 
gauge interactions. In order to systematically improve the result beyond the approximation 
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of the leading Yukawa terms, both effects from the gauge interactions and the momentum 
dependence should be taken into account. 

In the limit where the momentum dependence and the gauge interactions are neglected, 
the only missing contributions are the mixed two-loop Yukawa terms, 0{atab), and the 
0{al) corrections (and analogous contributions proportional to the Yukawa couplings of 
the other fermions and sfermions, which however are expected to give significantly smaller 
contributions than the third generation quarks and scalar quarks). As in the case of the 
0{abas) corrections, they can be numerically relevant only for large values of tan (3 and of 
the n parameter. Concerning the 0{atab) corrections, although their computation is not 
yet available, it is plausible to assume that the most relevant contributions are connected 
to the threshold effects in the bottom mass coming from chargino-stop loops JHl; which 
are actually already implemented in FeynHiggs as an additional option. An estimate of the 
0{al) corrections for a simplified choice of the MSSM parameters is presented in Ref. [T^ . 
According to that discussion, their effect can be at most comparable to that of the genuine 
0{abas) corrections, the latter being exemplified in Fig. |21 by the difference between the 
dashed and solid curves. 

To try to estimate, although in a very rough way, the importance of the various contri- 
butions we look at their relative size in the one- loop part. There, in the effective potential 
part, the effect of the 0{at) corrections typically amounts to an increase in rrih of 40-60 
GeV, depending on the choice of the MSSM parameters, whereas the corrections due to the 
electroweak (D-term) Higgs-squark interactions usually decrease rrih by less than 5 GeV j3] . 
Instead, the purely electroweak gauge corrections to rrih, namely those coming from Higgs, 
gauge boson and chargino or neutralino loops 0, are typically quite small at one-loop, and 
can reach at most 5 GeV in specific regions of the parameter space (namely for large values 
of n and M2). Concerning the effects induced by the dependence on the external momentum, 
as a general rule we expect them to be more relevant in the determination of the heaviest 
eigenvalue mn of the Higgs-boson mass matrix, and when M4 is larger than Mz- Indeed, 
only in this case the self-energies are evaluated at external momenta comparable to or larger 
than the masses circulating into the dominant loops. In addition, if is much larger than 
Mz, the relative importance of these corrections decreases, since the tree-level value of niH 
grows with M4. In fact, the effect of the one-loop momentum-dependent corrections on rrih 
amounts generally to less than 2 GeV. 

Assuming that the relative size of the two-loop contributions follows a pattern similar 
to the one-loop part, we estimate that the two-loop diagrams involving D-term interactions 
should induce a variation in rrih of at most 1-2 GeV, while we expect those with pure gauge 
electroweak interactions to contribute to rrih not very significantly, probably of the order 
of 1 GeV or less. Given the smallness of the one-loop contribution it seems quite unlikely 
that the effect of the momentum-dependent part of the 0{atas) corrections to rrih, which 
should be the largest among this type of two-loop contributions, could be larger than 1 GeV. 
As already said, the situation can, in principle, be different for the heavier Higgs-boson 
mass. However, the corrections to mn are relatively small in general and at one-loop level 
for most parts of the MSSM parameter space the momentum-dependent corrections are not 
particularly relevant. The momentum-dependent corrections turn out to be more relevant 
in processes where the H boson appears as an external particle, see Ref. [?7] . 

Another way of estimating the uncertainties of the kind discussed above is to investigate 
the renormalization scale dependence introduced via the DR definition of tan (3, fi, and the 
Higgs field renormalization constants. Varying the scale parameter between O.brrit and 2mt 
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gives rise to a shift in rrih of about ±1.5 GeV j2SI; in accordance with the estimates above. 

4.2 Estimate of the uncertainties from unknown three-loop cor- 
rections 

Even in the case that a complete two-loop computation of the MSSM Higgs masses is 
achieved, non-negligible uncertainties will remain, due to the effect of higher-order correc- 
tions. Although a three-loop computation of the Higgs masses is not available so far, it is 
possible to give at least a rough estimate for the size of these unknown contributions. 

A first estimate can be obtained by varying the renormalization scheme in which the 
parameters entering the two-loop part of the corrections are expressed. In fact, the resulting 
difference in the numerical results amounts formally to a three- loop effect. Since the 0{atas+ 
a^) corrections are particularly sensitive to the value of the top mass, we compare the 
predictions for rrih obtained using in the two-loop corrections either the top pole mass, 
^poie _ ]^';7^_2 GeV, or the SM running top mass m^, expressed in the MS renormalization 
scheme, i.e. 

pole 

m, - - l + 4„,(™,)/l-«,(m.)/2>r ' 

Inserting appropriate values for the SM running couplings and at we find mt = 
168.6 GeV. In Fig. we show the effect of changing the renormalization scheme for rrit 
in the two-loop part of the corrections. The relevant MSSM parameters are chosen as in 
Fig. m i.e. tan /? = 3 , Msusy = = /i = 1 TeV and nig = 800 GeV. The dot-dashed and 
dotted curves show the 0{atas) predictions for nih obtained using 171^°^'^ or frit, respectively, 
in the two-loop corrections. The solid and dashed curves, instead, show the corresponding 
0{atas + ttt ) predictions for nih. The difference in the two latter curves induced by the shift 
in rrit, which should give an indication of the size of the unknown three- loop corrections, is 
of the order of 1-1.5 GeV. However, as can be seen from the figure, the effect of the shift in 
rrit partially cancels between the 0{atas) and 0{a^) corrections and there is no guarantee 
that such a compensating effect will appear again in the three-loop corrections. 

An alternative way of estimating the typical size of the leading three-loop corrections 
makes use of the renormalization group approach. If all the supersymmetric particles (in- 
cluding the CP-odd Higgs boson A) have mass Msusy, and [3 = tt/2, the effective theory 
at scales below Msusy is just the SM, with the role of the Higgs doublet played by the 
doublet that gives mass to the up-type quarks. In this simplified case, it is easy to apply 
the techniques of Refs. [■19..20. in order to obtain the leading logarithmic corrections to rrih 
up to three loops (see also Ref. |SH1)- Considering for further simplification the case of zero 
stop mixing, we find: 

(Ami) 

(11) 

where t = log {M^^^y / m^) , rrit is defined in eq. (fTUj) . at and have to be interpreted as 
SM running quantities computed at the scale Q = mt, and the ellipses stand for higher 
loop contributions. It can be checked that, for Msusy = 1 TeV, the effect of the three-loop 
leading logarithmic terms amounts to an increase in ruh of the order of 1-1.5 GeV. If Msusy 
is pushed to larger values, the relative importance of the higher-order logarithmic corrections 
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Figure 5: ruh as a function of X^, using either mf°'^ or in the two- loop corrections. 
The relevant MSSM parameters are chosen as tan p — 3 , Msusy = Ma = h = 1 TeV and 
rrig — 800 GeV. For the different hues see text. 

obviously increases. In that case, it becomes necessary to resum the logarithmic corrections 
to all orders, by solving the appropriate renormalization group equations numerically. Since 
it is unlikely that a complete three-loop diagrammatic computation of the MSSM Higgs- 
boson masses will be available in the near future, it will probably be necessary to combine 
different approaches (e.g. diagrammatic, effective potential and renormalization group), in 
order to improve the accuracy of the theoretical predictions up to the level required to 
compare with the experimental results expected at the next generation of colliders. 

To summarize this discussion, the uncertainty in the prediction for the lightest CP-even 
Higgs boson arising from not yet calculatatcd three-loop and even higher-order corrections 
can conservatively be estimated to be 1-2 GeV. From the various missing two-loop corrections 
an uncertainty of less than 3 GeV is expected. However, it is extremely unhkely that all 
these effects would coherently sum up, with no partial compensation among them. Therefore 
we believe that a realistic estimate of the uncertainty from unknown higher-order corrections 
in the theoretical prediction for the lightest Higgs boson mass should not exceed 3 GeV. 
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5 Conclusions 



In this paper we have discussed the phenomenological impact of recently obtained results 
in the MSSM Higgs sector. The new corrections have now been implemented in the code 
FeynHiggs, which in this way provides the currently most precise evaluation of the masses 
and mixing angles of the CP-even MSSM Higgs sector. We have analyzed the effects of these 
new contributions by comparing the results obtained with FeynHiggsl.3, the newest version 
of the code that incorporates all these effects, with those derived with previous versions in 
which the two-loop 0{af) corrections were implemented via a renormalization group method 
and the sbottom corrections beyond one- loop were not taken into account. 

As a result of these improvements the lower limits on tan (3 that can be derived by com- 
bining the theoretical upper bound on rrih with its experimental lower bound, are significantly 
weakened. We have also shown that, if the top-quark mass is increased by one standard de- 
viation, the constraint on tan/? completely disappears if the other MSSM parameters are 
such that rrih takes its maximum value as a function of tan (3. 

Concerning the coupling of the lightest Higgs boson to the down-type fermions, that 
is related to the angle that diagonalizes the Higgs boson mass matrix including higher or- 
der corrections, we have shown that the area in the M4-tan /3-plane where its value gets 
suppressed by 70% or more with respect to the SM one is significantly modified. 

Finally, we have given an estimate of the uncertainties related to the various kinds of 
two- and three-loop contributions to the Higgs boson masses that are still unknown. Since 
it is extremely unlikely that all these effects would coherently sum up, the uncertainty in 
the theoretical prediction for the lightest MSSM Higgs boson mass from unknown two- and 
three-loop corrections should not exceed 3 GeV. 
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Note added: 

While finalizing this paper, Ref. jSH] appeared, containing the two-loop result for rrih obtained 
from the full two- loop effective potential of the MSSM [lU]. The contributions included in 
Ref. [Sn] that go beyond the ones discussed in this paper, namely momentum-independent 
contributions for non-vanishing gauge interactions and D-term Higgs-squark interactions, 
turn out to yield a shift in mh of 1 GeV or less in the numerical results given in Ref. PU] , 
thus confirming our estimate discussed in Sect. 14.11 
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